UNCLASSIFIED 


AD  NUMBER 


AD920489 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  only;  Test  and  Evaluation;  FEB 
1974.  Other  requests  shall  be  referred  to 
Air  Force  Armament  Lab,  Attn:  DLMA,  Eglin 
AFB,  FL  32542. 


AUTHORITY 


Air  Force  Armament  Lab  ltr  dtd  1  Apr  1976 


THIS  PAGE  IS  UNCLASSIFIED 


THIS  REPORT  HAS  BEEN  DELIMITED 
AND  CLEARED  FOR  PUBLIC  RELEASE 
UNDER  DOD  DIRECTIVE  5200,20  AND 
NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE, 

DISTRIBUTION  STATEMENT  A 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED, 


AD920489 


AFATL-TR-74-51 


Methods  For  Calculating  Magnus  Forces 
On  Slender  Bodies  Of  Revolution 


J.  E.  Milton 


Distribution  limited  to  U.  S.  Government  agencies  only; 
this  report  documents  test  ami  evaluation;  distribution 
limitation  applied  February  IP  K  Other  requests  for 
this  document  must  be  referred  to  the  Air  force  Armament 
Laboratory  (DLMA),  Eglin  Air  force  base,  Florida  32542. 


FOREWORD 


/  This  report  was  prepared  by  the  Eglin  Graduate  Center, 
College  of  Engineering,  University  of  Florida,  Eglin  Air  Force 
/Base,  Florida,  and  covers  work  performed  under  Contract  No. 

V  F08635-70-C-0065  with  the  Air  Force  Armament  Laboratory,  Eglin 
Air  Force  Base,  Florida.  The  project  monitor  for  the  Armament 
Laboratory  was  Dr.  George  B.  Findley  (DLMA).  The  work  described 
in  this  report  commenced  in  February  1972  and  was  completed  in 
December  1973. 


The  principal  contributor  was  Dr.  J.  E.  Milton. 

This  technical  report  has  been  I'eviewed  and  is  approved. 


RICHARD  M.  KELLER.,  Colonel,  USAF 
Chief,  Air-to-Surface  Modular 
Guided  Weapons  Division 


ii 


ABSTRACT 


Theories  for  predicting  the  Magnus  force  lor  slender, 
pointed  bodies  of  revolution  at  small  angle  of  attack  and 
small  spin  rates  are  reviewed.  A  semi-empirical  method  lor 
calculating  Cysp  is  discussed  and  shown  to  be  applicable  to 
a  variety  of  bodies  and  experimental  situations. 
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SECTION  I 


INTRODUCTION 


When  a  sphere,  cylinder,  or  other  body  of  revolution  which 
is  moving  through  a  fluid  is  rotated  such  that  the  axis  of 
rotation  is  at  an  angle  fro®  its  direction  of  translation,  an 
aerodynamic  force  is  produced.  This  force,  called  Magnus  force, 
is  in  addition  to  lift  or  drag  which  may  cccur  with  or  without 
spin.  The  direction  of  the  Magnus  force  is  perpendicular  to 
the  plane  formed  by  the  axis  of  rotation  and  the  translation 
direction.  The  magnitude  of  this  force  is  small  when  compared 
to  the  normal  force  which  may  be  developed  on  a  cylindrical 
body  typical  of  a  modem  missile.  It  is  found  to  be  less  than 
about  5  percent  of  this  normal  force.  However,  this  force  and 
its  associated  moment  about  the  body  center  of  gravity  may  be 
very  important  in  predicting  the  flight  path  of  spinning  pro¬ 
jectiles.  for  this  reason  it  is  of  considerable  interest  to 
aerodynamic ists  and  ballisticians.  Although  the  Magnus  effect 
has  been  observed  for  centuries  and  many  measurements  of  its 
magnitude  have  been  made,  its  origins  have  largely  not  been 
understood.  It  is  the  purpose  of  this  report  to  review  some 
of  the  theories  that  have  been  proposed  to  date  and  evaluate  them 
as  to  their  effectiveness  in  explaining  and  predicting  the  mag¬ 
nitude  of  this  force. 
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SECTION  II 


HISTORICAL  DISCUSSION 


According  to  Swanson^  the  first  record  of  the  observation 
of  the  drift  of  a  spinning  body  was  made  by  G.  T.  Walker  in 
1671,  the  body  was  a  sliced  tennis  ball.  L.ater  in  1730  B. 
Robins?  expressed  the  opinion  that  the  unexplained  dispersion 
of  cannon  balls  was  due  to  the  fact  that  they  were  spinning. 


In  order  to  try  to  put  the  ideas  on  a  more  firm  footing, 

G.  Magnus^  in  1852  performed  some  crude  experiments  both  in 
the  laboratory  and  in  the  field.  He  made  measurements  on  a 
spinning  cylinder  with  a  flow  normal  to  the  axis  of  the  cylinder 
and  found  that  unequal  pressures  were  developed  on  opposite 
sides  of  the  cylinder  due  to  the  combination  cf  the  normal  flow 
and  the  spin.  In  order  to  field  test  this  on  spheres,  he  fired 
musket  balls  which  had  the  center  of  mass  displaced  from  their 
geometric  center.  Then  by  placing  the  mass  center  to  the  right 
or  left  of  the  centerline  of  the  musket,  he  could  cause  a  spin 
to  the  right  or  left  to  develop.  He  found  that  in  this  way* 
he  could  cause  the  musket  ball  to  drift  either  to  the  left  or 
right  depending  upon  the  direction  of  the  spin.  This  demon¬ 
strated  that  the  earlier  ideas  had  been  correct,  and  the  phe¬ 
nomena  came  to  be  called  the  Magnus  effect. 

Lord  Raleigh4 ,  while  studying  the  flight  of  a  sliced  tennis 
all,  was  the  first  person  to  set  up  an  ideal  flow  representation 


of  the  situation.  He  demonstrated  that  if  one  superimposes  a 
uniform  potential  flew  on  a  circulation  a  force  in  the  proper 
direction  will  be  produced.  He  noted,  however,  that  his  model 
was  applicable  only  to  a  nor. -viscous  fluid.  His  theory  predict* 
no  £;€Cf<4ft&sm  5 or  th**  dove k-fMftt  of  the  needed  circulation  nor  « 
relationship  between  the  spin  rate  and  magnitude  of  the  force 
produced . 

These  results  wore  followed  by  experiments  by  La f ay  in 
Paris  between  19.0  and  This  work  was  dene  on  spinning 


Lif t .  'drag,  and 


cylinders  with  a  normal  flow  superimposed, 
pressure  distributions  wore  obtained-  for  two  different  rotors 
at  varying;  rotational  speeds,  Cne  significant  result  -of  this 
study  was  the  observation  that  at  lew  spin  *  a  negative 
Magnus  force  was  sometimes  observed.  This  »«ts  been  explained 
by  Krshn^c  and  is  associated  with  the  transition  of  the  boundary 
layer  from  laminar  to  turbulent  flow. 


was 


The  dilemma  of  hew  flu*'1  needed  ^es  row  lot  ion  could  be  developed 
solved  in  131*  when  !,.  Ft-andi  l7  shewed  that  it  could  result 


from  boundary  layer  vort icily  being  shed  by  separation.  The 
net  vorticity  of  the  shed  boundary  layer  was  equal  and  opposit 
to  the  resulting  circulation-  Recent  studies  by  Buford8,  Wixx  ■ 
and  Glauert*0  show  that  the  circulation  imparted  to  the  invis<  <1 
flow  is  always  less  than  that  o»  the  circulation  oi  the  me  . 
at  the  solid  surface. 

More  experimental  work  was  done  by  Reid**,  Thom**” and 
others  through  about  1935.  Since  that  time  many  tests  ». ave  been 
run  on  various  configurations  and  a  large  body  of  data  axist 
against  which  to  compare  any  proposed  theoretical  solution  to 
the  problem.  (See,  for  example,  references  18  through  23.) 

In  more  modern  times  a  number  of  people  have  tried  to  give 
a  solution  for  the  problem.  Examples  of  this  may  be  seen  in 
the  work  of  Martin2^,  Buford8,  Swanson*,  Kelley?8,  Kelley  and 
Thacker?6,  i’latou?7 ,  Iverson?8,  Rower29,  and  others. 

Generally,  the  theories  try  to  relate  the  Magnus  force  at 
high  angles  of  attack  to  the  problem  of  normal  flow  over  a 
spinning  cylinder.  In  most  cases  of  interest  ‘he  boundary  layer 
on  such  a  cylinder  will  have  separated,  and  the  mechanism  for 
the  product  ion  of  the  circulation  is  as  noted  before.  In  the 
case  of  typical  coni igurat ion  at  small  angles  a  different  mech¬ 
anism  is  needed. 


Martin’s?**  work  represents  a  step  in  the  direction  ot  trying 
to  make  predictions  by  the  use  el  modern  boundary  layer  theory. 

He  analysed  the  case  for  a  spinning  cylinder  In  incompressible 
laminar  flow  at  small  angle  of  attack.  He  obtained  a  solution 
tor  the  three-dimensional  boundary  layer  by  the  method  of  per¬ 
turbations.  With  this  he  determined  the  displacement  thickness 
of  the  boundary  layer  and  used  this  new  effective  shape  of  the 
body  t©  calculate,  by  slender  body  theory,  the  resultant  Magnus 
force.  This  same  approach  was  followed  by  Kelley-8  and  Kelley 
and  Thacker2®  who  extend  Mart  in ’ft*®  work  to  higher  npin  rates  and 
included  a  radial  pressure  gradient  and  sfcisv  friction  effects. 
These  theories  predict  a  linear  function  of  angle  'of  attack. 
Experiments  show  that  the  force  is  net  linear  with  angle  ei 
attack  even  at  lew  angles  and  shews  strong  non-linearity  at 
cross-flow  Reynolds  numbers  lower  than  where  one  would  expect 


to  i  ind  separation.  Turf  her,  exp»*t  ir-ent;.  by  Cover.  and  ftrkinr1 
arid  tHmtt3*  indicate  that  ter  acts -s§- inn 5 ng  N4 left  with  large  b/P 
and  at  low  angles  of  attack  the  tsoceta *  fetee  duo  to  croStsflew 
develops  with  distance  down  the  body  much  the  sure  as  the  flow 
dove  heps  about  an  -impulsively- traits  1  a  ted  circular  cylinder-. 
Kelley's32  work  on  the  prediction  of  the  norssal  force  on  a  sen- 
spinning  cylinder  at  an  angle  ©1  attack  used  this  approach. 


Platou^?  suggests  then  that  perhaps  the  Magnus  force  could  be 
treated  in  an  analogous  manner.  Platou  further  suggests  that 
the  normal  velocity  for  each  section  of  such  a  spinning  body 
be  taken  as  the  component  of  the  freestream  velocity  normal 
to  the  body  axis.  Further,  the  Magnus  force  on  each  section 
would  be  related,  not  to  the  steady  state  Magnus  forces  on  a 
cylinder,  but  rather  to  the  force  on  a  rotating  cylinder,  im¬ 
pulsively  rotated  from  rest  before  it  has  developed  the  steady- 
state  force. 

A  numerical  analysis  of  an  impulsively  rotated  cylinder 
immersed  in  a  uniform  freestream  was  performed  by  Thoman  and 
Szewczyk33.  From  this  work  the  cylinder  force  coefficient  can 
be  obtained  as  a  function  of  time.  Power  and  Iverson28,29  used 
Platou* s27  suggestion  and  the  results  of  Thoman  and  Szewczyk33 
to  predict  the  variation  of  the  Magnus  force  with  angle  of 
attack  and  spin.  The  results  predict  the  non-linearity  with 
angle  of  attack,  but  the  final  equation  contains  a  correlation 
constant  which  must  be  determined  from  experiment  and  which 
seems  to  take  on  a  large  range  of  values. 


'■  “1'  ‘  .  ,. .-$^J?Ar'J>'.~  .&*&r-:?t'i‘>.*!*  y‘.*l- £'■•■$  •'  ;. 


— .  -U.—'  -.-W-  c 


SECTION  III 
PROPOSED  THEORIES 


In  this  section  several  theories  will  be  discussed  that 
have  been  proposed  to  explain  the  Magnus  force  on  a  cylindri¬ 
cal  spinning  body  which  is  inclined  at  a  small  angle  a  from 
the  freestream.  The  boundary  layer  will  be  assumed  to  be  lam¬ 
inar,  and  no  provision  is  made  for  nose  shape.  Basically,  all 
of  these  theories  derive  from  the  work  of  Martin*’14  v,ho  was  the 
first  person  to  successfully  apply  modern  boundary  layer  theory 
to  this  problem.  Martin  reasoned  that  in  the  small  angle-of- 
attack  region  the  Magnus  force  would  not  be  due  to  circulation 
but  rather  to  a  distortion  of  the  effective  shape  of  the  body 
due  to  the  asymmetric  boundary  layer  which  would  develop  due 
to  t'  •»  spin.  If  this  shape  could  be  calculated,  potential  flow 
me+:  ; •  ic  could  be  used  to  determine  the  force.  Martin's  first 
paper  on  this  subject  was  in  a  cartesian  coordinate  system 
which  has  not  turned  out  to  be  the  simplest  one  to  use  for  this 
case.  In  addition,  there  were  some  errors  in  the  development. 
For  these  reasons  the  work  of  Kelley2'’ » 26  is  used  to  show  the 
basic  methods  for  this  solution  of  the  problem.  Kelley  intro¬ 
duced  the  cylindrical  coordinate  system  as  shown  in  Figure  1. 


The  Navier-Stokes  equations  for  steady  state  incompress¬ 
ible  flow  in  the  coordinate  system  shown  in  Figure  1  are: 
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where  u  and  v  are  the  velocities  in  the  x  and  y  directions  and 
w  is  the  velocity  in  the  9  direction  and  subscripts  denote  par¬ 
tial  differentiation. 


As  is  customary  in  Boundary 
tions  are  nondimensionaliscd  and 
they  are  reduced  to  the  boundary 
of  this  reduction,  see  reference 


Layer  Theory,  the  above  equa- 
by  an  order-of -magnitude  argument 
layer  equations.  (For  details 
2$.)  The  non-dimensional  boundary 
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Vj  £u«  1.  coordinate  System 
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,  «,  for  flow  over  a  cyx« 

layer  equation.  (5) 

may  be  written. 
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jj  =  u#  +  u,  +  u2  +  U3  +  u%  +  u5  +  u6 

J  =  v0  +  Vl  +  v2  +  v3  +  V,  +  v5  +  V6 

£  =  w3  w2  +  w3  +  w%  +  w5  +  w# 

where  u#  and  v#  are  non-dimensional  velocity  components  for 

zero  spin  and  zero  angle  of  attack.  The  subscript  1  refers 
to  terms  linear  in  a  or  w,  and  subscript  2  denotes  terms  of 
o.-der  a  ,  au,  to  ,  and  so  on  for  the  other  subscripts. 

To  zero  order  the  above  equations  reduce  to 

(rufl)x  +  (rv0)r  =  0  (b) 


u  u  ♦ 
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v  >1 
o  or 


-ft(ru 

Ur 


or^r 


(10) 


Equation  (4)  vanishes  since  for  a  =  w  =  0,  w  =  0  and  the  deriva¬ 
tive  with  respect  to  #  also  vanishes.  Using  numerical  integra¬ 
tion  techniques,  Seban  and  Dond^  solved  the  two  x'emaining 
equations  by  making  the  transformations 


The  solution  may  be  v.ritten 


(13) 


These  satisfy 
to  be  written 


Ux 


(nftl  -  f  -  Uc> 


(14) 


uation  (9)  identically  and  allow  equation  (10) 
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[ Cl+n5)fnr)3r|  ♦  ff„n  *  Uf£f, 


-  W  ■  0  • 


(15) 


This  can  be  reduced  to  a  system  of  ordinary  differential  equa¬ 
tions  by  assuming  f(n,0  to  be  a  polynomial  in  £: 


f(n,£)  =  f0(n)  +  fjCnH  +  f2(nU2  +  ••• 


(16) 


and  terminating  this  after  the  first  order  terms.  Using  this 
reduces  equation  (15)  to  three  third-ordei’  ordinary  differential 
equations  with  boundary  conditions  such  that 


u  =  v  =  0 


w  =  rw 


at  r 


uu-a  ) 


(18) 


UaSin^(l+~) 


for  large  r 


The  Navier-Stokes  equations  are  rewritten  six  times,  each 
time  considering  terms  of  one  higher  order  of  perturbation. 
Details  of  this  calculation  as  ivell  as  the  profiles  of  u,  v, 
and  w  are  presented  in  reference  26. 


Once  u,  v,  and  w  were  known 
how  far  a  streamline  at  the  cute 
would  be  moved  outward  from  the 
a  frictionless  fluid.  That  is  t 
method  for  doing  this  was  worked 
was  done  in  detail  by  Kelly  and 
tive  shape  induced  by  the  spin  s 


,  it  was  possible  to  calculate 
r  edge  of  the  boundary  layer 
position  it  would  have  had  in 
he  displacement  thickness.  The 
out  by  Dunn  and  Ke]ley3^  and 
Thacker*^.  This  gave  an  efl co¬ 
uch  as  shown  in  Figure  2. 


Slender  body  potential  flow  theory  was  applied  to  this, 
and  a  force  coefficient  was  obtained  in  terms  of  the  profile 
functions.  This  result  is  included  in  equation  (21). 

In  the  normal  thin  boundary  layer  theory  it  is  customary 
to  assume  that  there  is  no  pressure  gradient  normal  to  the  solid 
surface.  This  assumption  is  not  good  in  this  case,  equation 

(8)  shows  that  a  non-sero  value  for  M  will  be  obtained.  This 

*  3r 

equation  tc  order  £.  will  yield 

Li 


1  vr 

-  7 


(20) 
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Boundary  Layer 


Figure  2.  Effective  Shape 


Us  in*'  this  relationship  and  the  v  profile  functions,  Kelley 
and  Thacket’26  obtained  a  second  order  term  in  the  Magnus  force 
due  to  this  radial  pressure  gradient. 

As  is  well  known,  the  skiis  friction  shear  stress  at  a 
solid  boundary  such  as  this  problem  presents  is  given  by  Mew- 
ton's  relationship 

1  ■  *£L_  • 


if  this  is  integrated  over  an  area,  a  skin  friction  force  is 
obtained.  Taking  the  y  component,  this  can  be  shown  to  be 


2n 

Fy  =  /  r0xCos(J)d(J)  per  unit  length. 

y  0 

Since  the  gradient  lii  is  not  a  constant  around  the  circum- 

3r 

ference,  a  net  side  force  results.  Using  the  profile  functions, 
Kelley  and  Thacker28  obtained  a  third  order  contribution  to  the 
Magnus  force  from  this  mechanism. 

Combining  the  effect  of  displacement  thickness,  radial 
pressure  gradient,  and  skin  friction,  they  found 


Cy>pct 


[7.834  -  16.526 


+ 


.  (21) 


This  result  is  compared  with  a  number  of  experimental  measure¬ 
ments  in  Figures  3  through  23. 

As  could  be  expected,  due  to  the  assumptions  made  in  the 
derivation,  this  equation  does  not  give  good  results  above  very 
small  angles.  In  an  effort  to  rectify  this  situation,  Iver¬ 
son28  and  Power  and  Iverson2®  have  modified  this  scheme  to 
account  for  the  non-linearity  that  has  been  found. 


All  of  the  previous  work  has  assumed  that  no  circulation 
was  present  in  the  inviscid  outer  flow  since  there  was  no  separ¬ 
ation  of  the  boundary  layer.  It  is  known,  however,  from  work 
on  non-spinning  slender  bodies  that  at  fairly  small  angles  of 
attack,  about  six  degrees,  body  vortices  may  be  formed.  It 
seems  quite  plausible  then  to  assume  that  an  axial  distribution 
of  circulation  may  be  induced.  Power  and  Iverson2®  assumed  that 
this  circulation  could  be  approximated  by  assuming  that  the 
cross-flow  along  the  x-axis  of  the  spinning  cylinder  at  low 
angles  of  attack  is  similar  to  the  unsteady  flow  over  an  im¬ 
pulsively  rotated  cylinder  placed  perpendicular  to  the  f roe- 
stream.  The  boundary  conditions  at  the  outer  edge  of  the  bound¬ 
ary  layer  w‘ore  assumed  to  be: 


u  =  UC1-S.*) 

2 

"  =  u  Sinks'1) 

where  f 0  =  2*rev%.  and  k(x>  - 


and 


(22) 


(23) 
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where  r  is  the  circulation  that  actually  exists.  Note  that 
these  boundary  conditions  differ  from  the  ones  previously  used 
only  by  the  term  representin’'  the  shed  vorticity.  Values  of 
k(x)  were  found  by  using  the  results  of  Thoman  and  Szewc'/.yk*^. 
Power29  found  that  for  an  impulsively  rotated  cylinder  in  a 
crossflow  the  force  coefficient  was  well  represented  for  small 
time  by 


*  ■  ¥ 


(24) 


To  apply  this  to  the  spinning  cylinder  crossflow  let 
U_  =  USina 


(25) 


t  = 


UCosct 


(26) 


JU  = 

LL  bma 


(27) 


The  local  circulation  strength  is  given  by 

r  =  icjUj.d  . 


(28) 


Using  equations  (24)  and  (23)  one  can  write 


k(x)  o  x'  *  « 

r0  2* 


tana 


(29) 


which  induces  lor  small  a  to 


k(li)  *  OTI  Kc,s  C!  • 


C0  is  a  crossflow  correlation  constant  to  be  determined  from 
experiment  and  Rc  is  a  crossflow  Reynolds  number. 

With  these  new  conditions  Power^S  solved  the  boundary 
layer  equations  in  a  manner  similar  to  that  used  by  Kelley^. 
After  the  velocity  profiles  were  obtained,  contributions  to 
the  Magnus  force  by  displacement  thickness,  radial  pressure 
gradient,  skin  friction,  and  circulation  distribution  were 
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obtained.  The  final  result  is  given  by 


C  -  -  -ft- 

C7,P  * 


W 


60.97+( ,637R^75C0-357.37) 


1 

7E 


.227.01 

E 


V2 

vd; 


(30) 


Results  from  the  use  of  this  equation  are  compared  with  ex¬ 
perimental  test  results  in  Figures  17,  18,  and  22. 

Iverson?**  took  an  approach  to  the  problem  that  is  based 
on  the  same  basic  idea  as  the  work  described  above.  Instead 
of  going  through  a  boundary  layer  calculation  he  took  a  phe¬ 
nomenological  approach  and  tried  to  find  a  correlation  parameter 
directly  from  experimental  results.  He  used  equations  (25) 
through  (27)  repeated  below 


Up  =  USinct 
x 


t  = 


(25) 

(26) 


UCosa 


Sma 

and  equation  (24)  in  the  form 

Ci  =  fJ 

Putting  these  together  it  is  possible  to  obtain 


(27) 


,  s-1  1 

uct  (Vs^ 

Ml  f  l 

d  [Uc) 

M  ) 

Cl  = 


-1 

XT 


<p  tana  ( Ud }  S 


(31) 


This  is  a  local  side  force  coefficient  based  on  the  projected 
side-area  and  the  cross-flow  dynamic  pressure.  In  order  to  get 
it  for  the  complete  body,  simply  integrate  from  0  to  L  ami  change 
the  reference  area  to  wd*/4  and  the  dynamic  pressure  to  free- 
stream  to  obtain 


.p  s  kf(?)!  tano  sini/i,°(^) 


1 


(32) 


Iverson- 6  then  plotted  the  Magnus  force  coefficient  versus 
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2(lV  tana  Sin3^4a 

?ldJ 


.1 


for  a  number  of  experimental  cases.  He  was  thus  able  to  find 
k  which  according  to  the  theory  should  be  a  constant  good  for 
many  cases.  A  summary  of  results  based  on  this  equation  is 
given  in  the  Section  IV. 
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SECTION  IV 


RESULTS 


The  theories  outlined  in  the  previous  section  are  compared 
with  expei’imer.tal  results  of  a  number  of  investigators  working 
in  many  different  tunnels.  All  of  the  models  were  slender, 
pointed  bodies  of  revolution  with  no  boattailing.  Details  as 
to  Mach  number,  spin  rate,  and  such  are  given  on  the  figures. 

The  theory  that  has  been  commonly  used  to  predict  the 
Magnus  coefficient  for  slender  bodies  of  revolution  at  low 
angles  of  attack  is  that  of  Kelley2®.  The  result,  obtained 
before  the  work  with  Thacker28,  was  the  same  as  equation  (21) 
except  that  it  did  not  contain  the  higher  order  term  for  spin. 
The  result  is  given  by 


cy»pa 


62.672 


(33) 


It  is  found  that  the  coefficient  does  decrease  with  spin 
rate  but  not  at  the  rate  predicted  by  equation  (21).  Equation 
(33)  then  has  proven  to  be  more  accurate  than  the  higher  order 
one.  For  that  reason  the  results  presented  in  the  following 
figures  uses  equation  (33). 

The  results  of  the  work  by  Power  and  Iverson29  are  re¬ 
presented  by  equation  (30)  and  is  shown  for  comparison  on 
Figures  17,  18,  and  22,  The  value  of  C(J  that  was  used  is  shown 
on  each  figure. 

Equation  (32)  represents  the  semi-empirical  result  obtained 
by  Iverson28.  He  found,  on  looking  at  experimental  data,  that 
k  took  on  values  from  S  to  20.  Since  k  is  a  constant  of  pro¬ 
portionality,  the  Cytg  values  predicted  could  be  different  by 
as  much  as  a  factor  of  In  this  study  a  large  amount  of  ex¬ 
perimental  data  was  examined  in  order  to  find  a  reasonable 
value  of  k.  It  has  long  beer*  observed  that  the  Magnus  force 
increased  with  spin  rate  in  a  more  or  less  linear  fashion  for 
low  spin  rates.  At  higher  spin  rates  there  is  a  change  in  the 
slope  of  the  curve.  Since  equation  (32)  contains  no  provision 
for  accounting  for  this  non-linearity,  only  wind  tunnel  data 
in  this  spin  range  was  examined  in  this  attempt  to  evaluate  k. 
Based  on  these  restrictions,  a  value  of  k  =  10  was  selected. 

The  results  tor  this  are  shown  in  Figures  3  through  23. 
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SECTION  V 


DISCUSSION  AND  CONCLUSIONS 


As  is  obvious  from  Figures  3  to  23,  equation  (33)  gives 
reasonable  results  only  in  the  case  of  very  small  angles  of 
attack.  This  is  not  surprising  since  it  was  derived  for  a 
laminar  boundary  layer  in  incompressible  flow  with  no  separa¬ 
tion.  The  works  of  Martin?*4,  Kelley  and  Thacker^S >26  Were  an 
attempt  to  solve  the  problem  from  first  principles  and  were 
intended  to  be  the  forerunner  for  theoretical  work  on  more 
practical  cases.  Unfortunately  the  mathematics  involved  is  so 
complex  that  these  much  needed  improvements  have  yet  to  appear. 
Indeed  it  will  be  necessary  to  calculate  the  boundary  layer 
separation  lines  on  a  spinning  inclined  body  of  revolution  in 
order  to  get  the  solution  from  first  principles.  Some  idea 
of  the  magnitude  of  the  problem  may  be  obtained  by  realizing 
that  the  state  of  the  art  is  such  that  this  information  is 
only  beginning  to  be  available  for  non-spinning  bodies.  These 
results  are  still  open  to  serious  question  and,  in  fact,  the 
definition  of  what  constitutes  separation  of  a  three-dimensional 
boundary  layer  is  still  under  debaters. 

For  these  reasons  Power  and  Iverson?^  have  taken  a  semi- 
empirical  approach  based  on  suggestions??  that  the  problem 
could  be  treated  in  a  manner  similar  to  methods  used  to  pre¬ 
dict  the  non-linear  variation  of  forces  on  non-spinning  bodies. 

Power's  result,  equation  (30),  does  give  the  proper  shape 
of  the  Cy,p  versus  a  curves  as  seen  in  Figures  17,  18,  and  22. 

It  was  necessary  to  use  values  of  2.4  to  8  for  CQ  in  order  to 
make  the  equation  fit  the  data  in  just  these  three  cases. 

Neither  Power  nor  this  author  has  been  able  to  obtain  C0  in  a 
systematic  way  and  until  this  can  be  done  equation  (30)  will 
not  be  useful  in  predicting  the  Magnus  force. 


Equation  (32)  with  k  =  10  does  give  reasonable  answers 
for  most  of  the  cases  studied  in  this  report.  There  does  seem 
to  be  a  systematic  error  that  develops  in  the  transonic  region. 
This  is  not  surprising  since  the  prediction  of  any  aerodynamic 
characteristic  in  this  region  is  difficult.  Care  should  bo 
taken  in  the  application  of  these  results.  They  have  been 
demonstrated  to  be  valid  only  up  to  f?  =  0.4  and  for  straight 
afterbodies  with  no  fins  or  boattails.  A  variety  of  nose  con¬ 
figurations  are  represented  in  the  experimental  data:  cones, 
tangent-ogives,  Haaok-Hears ,  and  ethers.  At  these  low  angles 
ol  attack  the  Magnus  force  would  not  be  expected  to  depend 
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strongly  on  the  particular  shape  chosen^ ,  and  indeed  this 
appears  to  be  the  case  for  the  shapes  studied.  The  blunt  or 
hemispherical  nose  was  not  examined. 

It  is  the  opinion  of  this  author  that  the  non-linear 
variation  of  Cy  p  with  a  in  this  small  angle  region  is  assoc¬ 
iated  with  the  formation  of  the  so-called  body  vortex.  This 
begins  near  the  base  at  angles  of  attack  of  about  six  degrees 
and  as  a  increases  it  moves  forward  until  it  finally  reaches 
the  shoulder  of  the  nose^7.  if  this  were  true,  one  would  ex¬ 
pect  boattailing  to  have  a  strong  influence  on  the  behavior 
of  Cy  p.  This  is  indeed  the  case  and  has  been  previously 
reported  by  Platou^7  who  states  that  a  slight  rounding  of  the 
corner  of  the  base  caused  drastic  changes  in  the  Magnus  force, 
even  in  some  cases  changing  the  sign. 


Conclusions 


1.  Equation  (32)  with  k  =  10  will  yield  Cy>p  accurately 
enough  for  preliminary  design  calculations  for  pointed,  slender 
bodies  of  revolution  at  small  angles  of  attack,  with  p  <  0.4 
and  outside  the  transonic  region. 

2.  Much  work  remains  to  be  done  before  a  method  for  cal¬ 
culating  Cy<p  based  on  first  principles  is  devised. 

3.  In  the  small  angle-of -attack  region  the  Magnus  force 
does  not  depend  strongly  on  the  particular  pointed  nose  that 
is  chosen. 

4.  More  work  should  be  done  in  order  to  extend  the  appli¬ 
cability  of  equation  (32)  into  the  region  where  Cy>p  is  a  non¬ 
linear  function  of  p  and  if  possible  into  the  transonic  region. 
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Theories  for  predicting  the  Magnus  force  for  slender,  pointed 
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rates  arc  reviewed.  A  semi-empirical  method  for  calculating 
Cv,p  is  discussed  and  shown  to  be  applicable  to  a  variety  of 
bodies  and  experimental  situations. 
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